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Abstract 

All high-temperature phases of the known N = 4 superstrings in five dimensions can be 
described by the universal thermal potential of an effective four-dimensional supergravity. 
This theory, in addition to three moduli s, t, u, contains non-trivial winding modes that 
become massless in certain regions of the thermal moduli space, triggering the instabilities 
at the Hagedorn temperature. In this context, we look for exact domain wall solutions 
of first order BPS equations. These solutions preserve half of the supersymmetries, in 
contrast to the usual finite-temperature weak-coupling approximation, and as such may 
constitute a new phase of finite-temperature superstrings. We present exact solutions for 
the type-IIA and type-IIB theories and for a self-dual hybrid type-II theory. While for 
the heterotic case the general solution cannot be given in closed form, we still present a 
complete picture and a detailed analysis of the behaviour around the weak and strong 
coupling limits and around certain critical points. In all cases these BPS solutions have 
no instabilities at any temperature. Finally, we address the physical meaning of the 
resulting geometries within the contexts of supergravity and string theory. 



1 Introduction 



In the conventional description of string theory, where one has only very little information 
about the dynamics of all possible string states, the main framework is provided by an 
effective supergravity for the massless fields which is obtained by integrating out the 
massive modes. As a result, many stringy effects that are attributed to the massive 
modes cannot be addressed in a systematic way; we can only appreciate their relevance 
in certain regions of the moduli space where massive states can become massless. In those 
cases, the conventional supergravity approach has to be enlarged to include the massless 
as well as all the relevant would-be massless states on an equal footing, for otherwise the 
effective theory will break down due to the appearance of singularities. Put differently, 
the possibility to have extra massless states in string theory signals the limitations of 
the perturbative field theory description of string dynamics in all corners of the moduli 
space. Conversely, including the would-be massless fields into the effective field theory 
can teach us important non-perturbative lessons about string theory. 

There are specific interesting problems which can be addressed systematically, based 
on supersymmetry, by isolating the dynamics of a few relevant modes that can become 
massless. An important example of this kind is provided by the conifold singularity in 
field theory and its string-theoretical resolution Jj], ^ . In particular, in Calabi-Yau com- 
pactifications which admit a non-trivial 3-cycle, the moduli space develops a singularity 
that invalidates the applicability of the conventional low-energy effective theory when 
the cycle shrinks to zero size. Incorporating a D-brane that wraps the 3-cycle, with mass 
proportional to its period, resolves the problem of conifold singularities as the enlarged 
theory is appropriate for describing cycles of all sizes. 

Superstrings at finite temperature provide another interesting example that has been 
studied on and off for quite some time, but admittedly has not been fully explored 
yet. In particular, it is known that string theory exhibits an exponential growth of the 
number of states at high energy, which gives rise to a limiting temperature, known as 



the Hagedorn temperature Th |3]|— 1|11|| . From the world-sheet point of view, where one 
uses a periodic Gaussian model to describe the propagation of strings on tori, there is 
a Kosterlitz-Thouless phase transition at a critical radius where vortices can condense 



I~2| , |X3 1 , thus leading to a limiting temperature for string thermodynamics (see also 



TT| ) associated with a phase transition. It has been further realized that there are string 
states with non-trivial winding number that can become light close to this temperature, 
and then turn tachyonic beyond it, thus signaling thermal instabilities of string theory 
at very high temperatures ||14|| . Some implications of this effect have been examined 
in the context of strings in the very early universe, where one tries to find ways that 
avoid the initial singularity, and at the same time explain why the dimensionality of 
the physical space-time is four ]T3| (see also and |I7| for a more recent exposition 
that takes into account D-branes). In an recent development, the contribution of the 
relevant winding modes was explicitly described in terms of an effective four-dimensional 
supergravity using a universal thermal potential that incorporates all phases of N = 4 
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superstrings, taking into account the thermal deformation of the BPS mass formula in 
N = 4 supersymmetric theories JTB], [L9|- However, no solutions have been found so far 



with a concrete physical interpretation and the ability to resolve the problems of quantum 
cosmology or the subtle issues raised at the end point of the black-hole evaporation; these 
are situations where the temperature grows to infinity and one encounters singularities 
in the context of any perturbative field theory. 

The present paper grew out of the attempt to construct explicit solutions of the 
effective theory of supergravity which was proposed to include the dynamics of the lightest 
relevant winding modes that could become massless at the Hagedorn temperature. We 
construct for the first time families of non-trivial solutions with varying winding fields 
in the type-II and the heterotic sectors of string theory. We focus attention primarily 
on domain wall solutions, one of the reasons being that they are technically the simplest 
ones to consider; all the fields in the bosonic sector of the theory, namely the metric and 
the collection of scalar fields, are taken to depend on a single spatial coordinate only. 
Ultimately, of course, it is necessary to search for other types of solutions, for instance 
solutions with spherical symmetry, in order to discuss models appropriate for quantum 
cosmology or black-hole physics. We hope to report on such generalized solutions of the 
effective supergravity elsewhere in the future. 

There are also some other compelling reasons for being interested in the existence 
and the explicit construction of domain wall solutions of the effective supergravity that 
describes strings at finite temperature. The first is provided by the important role of 
domain walls, as space-time defects, in a variety of cosmological applications, though we 
do not focus here on these types of physical problems. The second is provided by the fact 
that the domain walls, by their nature, break the four-dimensional Poincare invariance to 
an effective three-dimensional invariance. Then, one can have a working framework for 
studying the issue of supersymmetry breaking by non-perturbative effects, as it was done 
in fl8], pLQfl , taking into account the peculiarities of supersymmetry in three dimensions, 
namely that even if local supersymmetry is unbroken the massive multiplets will not 
necessarily exhibit mass degeneracy [20]. 

The domain wall ansatz allows us to consider solutions of certain first order differ- 
ential equations derived from a prepotential. As we explicitly show, these solutions by 
construction preserve half of the supersymmetries, and in this sense they are BPS. The 
full set of equations allows consistent truncations to subsectors which correspond to the 
heterotic string, type-IIA or type-IIB, or the type-II string at the self-dual radius (which 
we call the hybrid type-II). For all type-II cases we can find the exact general solutions in 
closed form. For the heterotic sector, however, no solution in terms of known functions 
is available, and we instead give a detailed analysis of the solutions around the weak and 
the various strong coupling limits, as well as around certain critical points, which still 
allows us to obtain a reasonably complete overall picture. 

It should be stressed that in the standard finite-temperature treatment of perturbative 
superstrings all supersymmetries are broken due to the different boundary conditions 
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on bosons and fermions one has to impose in the periodic imaginary time direction. 
For this reason, it is not surprising that some modes become tachyonic beyond certain 
temperatures. On the contrary, the solutions we find within the domain wall ansatz of 
the effective supergravity are non-pertubative, containing regions of strong coupling, and 
preserve half of the supersymmetries. As such they are expected to be stable solutions, 
and indeed we show that, although the temperature can be arbitrarily high, no tachyonic 
instabilities ever develop. Even though we are only considering solutions of an effective 
supergravity, they may well point to a new finite-temperature phase of superstrings which 
is BPS and has no thermal instabilities, i.e. no Hagedorn temperature. 

To further probe the physical meaning of these domain wall solutions we first study the 
propagation of a quantum test particle in the corresponding supergravity backgrounds. 
This allows to discriminate between wave-regular and non-regular geometries, thus nar- 
rowing down the physically meaningful set of solutions. The analogous question could 
be addressed within string theory. The consistency of string thermodynamics seems to 
require the compactness of all spatial directions since otherwise one would encounter 
problems with negative specific heat (see, for instance, [15] for a delicate analysis, with- 
out making reference to the geometrical details of the space, and references therein). It 
is clear from our explicit solutions in the type-II cases, and also true for the heterotic 
case (with the possible exception of certain solutions that are never weakly coupled), 
that none of them is periodic in the spatial dimension that defines the domain walls. 
However, this does not mean that our solutions of the effective supergravity theory are 
insufficient to extrapolate to string theory. Indeed, the compactness criterion is estab- 
lished within the micro-canonical treatment of string thermodynamics, and this assumes 
that the theory remains weakly coupled throughout all space. Furthermore, it concerns 
the phase where supersymmetry is completely broken by the finite temperature. On the 
other hand, almost all our solutions always contain a region of strong coupling. Most 
important, our solutions are BPS, preserving half of the supersymmetries. These differ- 
ences are enough to cast serious doubt on whether the arguments about compactness of 
all spatial dimensions should apply. Furthermore, as already noted, the fields correspond- 
ing to the winding modes of our solutions never become massless even if the temperature 
modulus approaches the would-be Hagedorn temperature, which is yet another way to 
see that these solutions belong to a different phase. 

This paper is organized as follows: In section 2, we review the main aspects of the 
effective theory of supergravity that was proposed to describe the thermal phases of all 
N = 4 superstrings and present the form of the universal thermal potential for the six 
scalar fields that correspond to the s, t, u moduli and the winding states zi, Zi, z 3 that can 
become massless in the heterotic, type-IIA and type-IIB sectors. In section 3, we consider 
the domain wall ansatz for the metric and all the relevant scalar fields. We first discuss 
quite generally when a D = 4, N = 1 supergravity theory admits domain wall solutions 
that are derived from first order differential equations. Indeed, this follows from some 
simple reality assumptions about the scalar fields and the Kahler potential. We show that 
these solutions automatically preserve half (or all) of the supersymmetries and hence are 
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BPS. Then we specialise to the effective thermal supergravity and explicitly derive the 
coupled system of six first order non-linear differential equations for them. Consistent 
truncations of this system lead to various type-II or heterotic sectors. In section 4, this 
system is analysed explicitly in the type-IIA, type-IIB as well as in the hybrid type-II 
sectors where the general solutions can be derived in closed form. In section 5, we study 
the equations for the heterotic sector which cannot be solved analytically, apart from a 
very special solution. We are able to study the behaviour of the solutions in the vicinity 
of the weak and the strong coupling regions, as well as around certain critical points, 
and present a fairly complete general picture. The construction of domain wall solutions 
for the whole system of six equations is beyond the scope of the present paper and can 
probably only be done numerically. In section 6, we study the role of boundary conditions 
in selecting physical solutions in the effective theory of supergravity and comment on 
string thermodynamics. Finally, in section 7, we present our conclusions and outline 
some directions for future work. In an appendix, we address the issue of having periodic 
solutions using criteria from the general theory of dynamical systems. This may turn out 
to be useful when studying the complete system of six equations. 



2 Thermal potential of effective supergravity 

In this section we present the essential ingredients for constructing an effective N = 1 
supergravity in four dimensions^ that describes the thermal phases of all N = 4 super- 
strings in a universal way, following earlier work [|18|, [2lJ. The construction takes 
into account the dynamics of the would-be tachyonic winding modes of the N = 4 su- 
perstrings which are responsible for inducing a phase transition at high temperatures 
in string thermodynamics. As such, the effective N = 1 supergravity theory provides 
a systematic framework for quantifying arguments about the phase transition occuring 
at the Hagedorn temperature due to the special form of the effective potential of the 
relevant winding modes u, which includes a trilinear coupling auui with the modulus a 



describing the temperature 14 



The starting point of the construction is provided by five- dimensional N = 4 theories 
that are effectively four-dimensional at finite temperature. The crucial observation is that 
ci-dimensional superstrings at finite temperature look like (d— l)-dimensional strings with 
spontaneously broken supersymmetry. Concretely, putting strings at finite temperature 
T amounts to compactifying the (Euclidean) time on a circle of inverse radius R" 1 = 2nT . 
As the Euclidean time direction is compact, one imposes boundary conditions to take 
particle statistics into account; modular invariance of the thermal partition function 
dictates then specific phase factors in the related GSO projection fT4|. Technically, this 



procedure is equivalent to a Scherk-Schwarz compactification from d = 5 to d — 1 = 4 



1 We always count four-dimensional supersymmetries, i.e., N = 1 supersymmetry has four super- 
charges. 
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with a well-defined gauging associated to the temperature modification of the effective 
theory of gauged supergravity |^2, 



1 8 1 . This gauging of N = 4 spontaneously breaks 



super symmetry 



2.1 Thermal dyonic modes 

We restrict ourselves to the study of N = 4 strings, because in this case the supersymme- 
try algebra and its central extensions suffice to determine the masses of all BPS states. 
Recall that the states which become tachyonic above a certain Hagedorn temperature 
Tfi are 1/2-BPS states, at least in dimensions where there are no other BPS states with 
smaller fractions of supersymmetry. Then, in N = 4 strings, one can identify completely 
all the perturbative and the non-perturbative BPS states that can induce thermal insta- 
bilities and construct from first principles an exact effective supergravity for these states, 



as it was done in refs. JL81 [19|] . It is useful to observe that the odd dyonic modes exhibit 
the same finite-temperature behaviour as the odd winding string states. This follows 
from the action of dualities on the N = A strings, relating the dyonic modes to perturba- 
tive winding states of dual strings. These windings in turn induce thermal instabilities 
in the dual strings. Describing the complete set of thermal phases of N = 4 theories 
requires then both perturbative and non-perturbative states from all string viewpoints. 

We define the context of the present work, following |18|, [19] , by first considering string 



theories in six dimensions obtained by compactification on T 4 (heterotic string) or K3 
(type-II strings). We further compactify one dimension S l , with radius Rq, and so the re- 
sulting five-dimensional theory exhibits T-duality between type-IIA and type-IIB (under 
Rq — > a' hct /Rs) and S-dualities between heterotic and type-II strings. Moreover, since we 
are putting strings at finite temperature T, the time dimension is taken Euclidean and 
compactified on S 1 with radius R = l/(2irT) with twisted fermionic boundary conditions 
that account for the thermal effects. 

In order to obtain the resulting four- dimensional thermal mass formula and examine 
which states can become tachyonic above a certain temperature, thus inducing instabil- 
ities, we first consider the usual BPS mass formula in N = 4 supersymmetric theories 



24], £5], ^7j] written in heterotic variables: 
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a' het tu 



m + ntu + iiwbU + n't) + is[m + ntu — i(m'u + n't)} . (2.1) 



Here s, t and u are defined in terms of the compactification radii Rq, R and the heterotic 
string coupling ghet as follows, 



RRq r 

/ ' U = ~D 

9hct "het -^6 



, t = -^ } u= — , (2.2) 



supplemented by the relation a' het = 4s, in units where the four-dimensional gravitational 
coupling k has been normalized to y/2. These expressions will help us later to under- 
stand the physical meaning of the various domain wall solutions in terms of the three 
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physical parameters in the problem R 6 , R and phet- The integers m, n, mf, n' are the 
electric momentum and winding quantum numbers associated to the compactification 
on the 2-torus with radii R§ and R. The tilded integers are the corresponding magnetic 
non-perturbative counterparts. Then, string dualities are simply described by the inter- 
changes s <r-* t (heterotic-IIA), s <-> u (heterotic-IIB) and t *-> u (IIA-IIB), which leave 
invariant the mass formula and the temperature radius R 2 = 4stu which is common to 
all three string theories when measured in Planck units. 

As discussed in refs. [JTR IS], the thermal deformation of the mass formula ( J2.1[ ) is 
very simple: the momentum quantum number m is replaced by m + Q' +n/2, Q' being the 
helicity charge. It also reverses the GSO projection, to account for the modified boundary 
conditions in the temperature deformation of the theory. For our present purposes, we 
may restrict our attention to only the states having m' = n' = m = n = 0, since the 
states that can become tachyonic first as the temperature increases are contained within 



this subset ||19|| . The thermal spectrum of these light states is then given by 

M% = (m + Q' + ^kp^j + k T Ptq ^ r R^j - 2T W 8\ k \ A 8 Q/}0 , (2.3) 

where we have defined the integers k, p, q and r by the condition (n, mf, h') = k(p, q, r), 
p, q and r being relatively prime numbers, and T p ^ r denotes the effective string tension 

^ v q t , . 

T w = -P- + + — , (2.4) 

a het a llA a llB 

which can be written in terms of s, t and u using the identifications 

"Let = 4 s , a'uA = 4t , a 'im = 4u , (2.5) 

again in the Planck units with normalization k = The mass formula (|2.3|) is an 
extension of the known perturbative finite-temperature mass formulas with the correct 
duality and zero temperature behaviour. In particular, the last term is generated in 
perturbative strings by the GSO projection related to the breaking of supersymmetry. 
The various constraints imposed among the quantum numbers require that p, q, r are all 
positive and that mk > — 1. 

Analysing the four-dimensional mass formula for Ai?, we note that tachyons can 
appear when Q' — and \k\ — 1 for certain values of the four quantum numbers m, 
n, m', h' . The critical temperature (and hence the critical radius R) can be found by 
locating the zeros of A4^. The first tachyons and the critical temperatures are: 



(m, n, rh' n') 



±(-1,1,0,0) with R={V2 + 1)\I^ , (2.6) 



(m,n,m,n) = ±(0,0,1,0) with R = j2a' IIA , (2.7) 



(m,n,m,n) = ±(0,0,0,1) with R = y2a' nB . (2.8) 

Since the winding numbers in the heterotic, the type-IIA and the type-IIB strings are 
respectively n, mf and n', these three pairs of states have winding numbers ±1, as ex- 
pected on general grounds, in their respective perturbative superstring theory. There 
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are also two other series of states occuring for Q' — and \k\ = 1, with quantum num- 
bers m = — 1 and (p, q, r) arbitrary, which can become tachyonic when p = 1 or p = 2; 
however, in either series the critical temperatures that result from M.\ = are higher 
than at least one of the perturbative Hagedorn temperatures of the heterotic, type-IIA 
or type-IIB strings given above. 

This analysis of the mass formula suggests the field content to be used in the effective 
field theory description. We certainly need the three moduli s, t and u and the six scalar 
fields (actually three conjugate pairs of fields) able to generate the thermal transitions. 
These states will be embedded in the generic scalar manifold of N = 4 supergravity, 

Sl(2,R) SO(Q,m) 



U(l) 50(6) x SO(m) ' 

where m is the number of vector multiplets. To study the thermal phase structure 
however, it is consistent and sufficient to truncate this N = 4 theory to N = 1 and to 
retain only chiral N = 1 multiplets describing the relevant moduli and winding states. 
After truncation, the effective theory includes three complex moduli scalars S, T and U, 
with scalar manifold 

( SL(2,R) \ ( SL(2,R) \ f SL(2,R) \ 

{-mTTJs {~uw~)r {~m~)u • (2 ' 9) 



and six complex windings Zj[, A = 1, 2, 3, living on 

( SO(2,3) \ ( SO(2,3) 
\S0{2) x S0(3)J z+ \S0(2) x 50(3) , 



(2.10) 



The resulting scalar manifold is a Kahler manifold for chiral multiplets coupled to N = 1 
supergravity that arises by two successive Z 2 projections applied on the N = 4 scalar 
manifold. 



2.2 The effective supergravity 

The construction of the effective N = 1 supergravity Lagrangian proceeds as follows 
li~9| . A generic four- dimensional N = 1 supergravity theory is characterized by a Kahler 



potential K and a holomorphic superpotential W. The bosonic sector of the theory is 
given by the following Lagrangian density: 

e~ x C = \R- ^Kjjid^id^) - V($, $) , (2.11) 

where Kjj = d 2 K/d§ I d§ J are the components of the Kahler metric for the collection 
of complex scalar fields {$ 7 } present in chiral multiplets. The scalar potential V($, $) 
assumes the special form 



V 



- e K (K 1 j W ;I W. j - 3WW) , (2.12) 
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using the notation of covariant derivatives 

dW OK 

w -' = w + w w - < 2 ' 13 > 

The form of the Kahler potential K follows from the constraints defining the N = 4 
scalar manifold, truncated to N — 1. One finds, in particular, 

K = -ln(S + S)~ ln(T + T) — \n(U + U)- lnF+ - \nY_ , (2.14) 

with 

Y ± = l- 2Z±Z ± A + {ZiZ±)(Z%Z%) , (2.15) 

where summation over repeated indices A or B = 1,2,3 is implicitly assumed. The 
expression of the superpotential W follows from the gauging applied to the N = 4 
theory. This gauging encodes the Scherk-Schwarz supersymmetry breaking mechanism 
generated by the introduction of a finite temperature |HJ and yields 

1 



w = 2V2 



-yl - Z\ZX){\ - Z B Z B ) + (TU - l)Z+Zr + SUZtZ 2 + STZtZ^ 



(2.16) 

The scalar potential of the theory defined by K and W is complicated, but focusing 
on the real directions defined by ImZ^ = ImS* = ImT = ImU = Re (Z\ — Z~^) = 
leads to important simplifications. Furthermore, a study of the mass spectrum in the 
low-temperature limit shows that these directions correspond precisely to the possible 
phase transitions occuring in the theory. We are then led to consider only the relevant 
would-be tachyonic states using the field variables 

s = ReS, t = ReT, u = ReU, z A = Re Z\ = Re Z A . (2.17) 

Let us introduce for later convenience the quantities 

.3 



A=l 1 x 



Then, in this case, the components of the Kahler metric in the field space become 

Kss = h> Ktt = ^> KuD = ir^ Ka±s± = o^? 61 > (2 - 19) 

i.e., the metric becomes diagonal in the directions Z A ; in writing the kinetic terms for 
the fields za we have to take K A b = 45^/(1 — x 2 ) 2 , as there is a factor of 2 picked up by 
Z\ = Zj in this case. 

In order to present the Lagrangian density for all remaining six scalar fields, we find 
it helpful to trade s, t and u for (pi, (f> 2 and </> 3 as follows: 

s = e" 201 , t = e~ 2 ^ 2 , u = e" 2 ^ 3 . (2.20) 



S 



Note that fa — > +00 corresponds to the strong coupling limit s — * 0, whereas the weak 
coupling limit is attained for (f>i — ► —00. With these definitions, we obtain a simplified 
form of the effective supergravity with bosonic Lagrangian density 

^ = \ R -\ (g(^) 2 + E jr^r^ 2 ) - v ' ( 2 - 21 ) 

where the scalar potential is now given by the sum 

V = V 1 + V 2 + V 3 . (2.22) 
Each term Va is a function of the moduli 0j and a polynomial in if^. Explicitly, we have 

Vi = ^e 2 * 1 if 1 2 (cosh(20 2 + 203)(4if 1 2 + l)-3) , (2.23) 
V 2 = ^e 2 * 2 if 2 2 (e^ 1 " 2 ^(4F 2 2 + 1) - 4) , (2.24) 
V 3 = -e 2 ^Hi (e- 2 ^- 2< ^(4# 2 + 1) -4) . (2.25) 
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Our normalizations have been chosen so that the kinetic terms of the fields fa assume 
their canonical form, with an overall factor 1/2 as coefficient, since the normalization 
k = v^2 has been used. Then, the scalar potential V provides the universal thermal 
effective potential that describes all possible high-temperature instabilities of N = 4 
superstrings. 

At this point, it is very useful to observe that the scalar potential ( |2.22j ) for the six 
real fields s, t, u and za can also be cast into the special form 

This property only depends on the structure of the Kahler potential K. The 'prepotential' 
W given by 



(2.27) 



W = I e ^ 1+ ^+* 3 - 2e^sinh(0 2 + <p 3 )Hf + e'* 1 (e^'^Hl + e^'^Hl 

= 7k (I + ( tu - l ) R l + suH 2 + ■ 

The second expression indicates that the prepotential is actually given by 

W = e K/2 W I 

real directions i 

(2.28) 

as can be inferred from the general form of the supergravity potential and from our 
specific Kahler potential^ fl2.14j) and superpotential (|2.16|) . As we will see in section 
3.2, this is not an accident. The prepotential will be of particular importance when 
writing the system of first order differential equations for the domain wall configurations 



2 In particular, exp(— K/2) = 2\j2stu (1 — x 2 ) 2 upon restriction to the real directions of the theory. 
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of the theory. It should be emphasized that given the potential (|2.22 ), we may view 



( p.26| ) as a differential equation for the prepotential W. Hence, in general, there can be 
solutions other than (|2.27|) . However, the latter choice is the one that leads to one-half 
super symmetric solutions, as we will see in section 3.3. In general, another choice for W 
satisfying ( |2.26| ) would break supersymmetry completely. 



2.3 Thermal phases of N = 4 strings 

The thermal phase structure of N = 4 superstrings can be deduced from the form of the 
potential V. First of all note that the low-temperature phase, i.e. large R, or large stu, 
or large and negative 01 + 02 + 03, is common to all strings. It is characterized by the 
condition Hi = H 2 = H 3 = 0, in which case V\ = V 2 = V 3 = 0, and so the potential 
vanishes for all values of the moduli fields s, t and u. Since the four-dimensional couplings 
of the three strings are f\ 

V2 \/2 r RR G V2 ^R , 

s = -V- , t = = , u = \- = V2— , (2.29) 

9het 9llA a het 9llB K 6 

we see that this phase exists in the perturbative regime of each string theory. Notice also 
that the configuration Ha = and 0j constant (but arbitrary) is an exact solution of the 
theory, with all super ■symmetries broken by the thermal deformation. The reason is that 
in this phase the supersymmetry transformation is proportional to 

e K/2 W (2 30) 

which does not vanish in the moduli directions, except of course in the limit stu — > oo 
(zero-temperature limit). 

There are three high-temperature phases generated by each of the three contributions 
to the scalar potential Q2.22 ), and which are characterized by a non-zero value of a winding 



state za- Which phase is selected depends on the values of the moduli i; corresponding 
to the coupling, the temperature radius R and the sixth-dimension radius _R 6 in each 
string theory (with string units a' het , ot' llA or a' I1B ). None of them is a solution to the 



theory with constant scalars. The detailed analysis of the thermal phase structure ||19 | 
leads to the following classification. 

Het erotic high-T phase: 

In heterotic-string units, the temperature modulus is tu = l/(27r 2 a[ iet T 2 ), while it is 
stu = 1/(27tkT) 2 in Planck units. A heterotic tachyon Hi appears first if 



l) 2 < tu < (V2 + l) 2 , su > 4 , st > 4 , (2.31) 



3 Notice that there is a factor of \/2 that rescales the right hand sides of (2.2E) as compared with the 



corresponding equation for zero temperature (2.2). This is due to the new spin structure introduced by 



temperature. However, eq. (|2.5| ) for the three string scales remains the same. 
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where the last two condition guarantee the absence of type-II tachyons. Then, this implies 
for the four-dimensional heterotic coupling that 

V2S- 1 = g 2 lct < g 2 Tit = . (2.32) 

The upper and lower critical temperatures in ( |2.31| ) are related by T-duality in the 
temperature radius. Within the interval (|2.31|) , the theory has a (non-tachyonic) solution 
for Hi = 1/2, tu = 1. And, since the scalar potential reduces to 

V = ~e 2(t>1 , (2.33) 

a linear dilaton background is a solution. This background breaks one half of the super- 
symmetries [PJ. This particular solution will be rederived as a special case of a more 
general class of domain wall backgrounds that break half of the supersymmetry in the 
heterotic sector, in section 5. 

In the strong coupling regime of the heterotic theory, where s — > 0, there are only 
type-II instabilities; the high-temperature heterotic phase cannot be reached for any value 
of the radius R 6 if g^ et > g 2 T if Put differently, the high-temperature heterotic phase can 
only exist in the perturbative heterotic regime, which by the heterotic-type-II duality 
corresponds to the non-perturbative regime of the type-II string theory. 

Type-IIA/B high-T phases: 

The relevant temperature modulus is either su = l/(2n 2 a' llA T 2 ) for type-IIA or st = 
1 / (2n 2 a' im T 2 ) for type-IIB strings. The high-temperature phase for type-II strings is 
then defined by su < 4 (for IIA) or st < 4 (for IIB). Thus, type-II instabilities arise in 
the strong coupling region of the heterotic string. The analysis of the problem shows that 
the value Rq/^J oe' het determines whether a type-IIA or a IIB tachyon comes first as the 
temperature increases. We have in particular the following two high-temperature type-II 
phases: 

R, < , (2.34) 

R 6 > (2.35) 

Note that the type-II high-temperature phases differ from the heterotic one by the 
absence of an exact solution with constant windings and/or moduli. In these phases, the 
theory depends on a non-zero, field-dependent winding and two non-trivial moduli, t and 
su in the IIA case, u and st for IIB strings. This concludes our general review of the 
thermal aspects of iV = 4 strings at finite temperature. 



IIA 



IIB 



2ttT > 



2ttT > 



2 fl'het ^6 

1 Rn 



2g 



het a het 



9het > 



9het > 



V2 + 1 

V2 + 1 
2V2 
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3 BPS domain wall solutions 



We are now in the position to look for solutions of the effective supergravity that describe 
the thermal phases of iV = 4 superstrings. We will use the domain wall ansatz for two 
reasons. First, this probably is the simplest possible class of solutions one may construct 
in supergravity theories whose bosonic sector is described by gravity coupled to a selection 
of scalar fields with non-trivial potential terms, as the universal thermal effective potential 
we have at our disposal. Second, with all fields depending only on a single variable one can 
look for solutions that satisfy more stringent, but easier to solve, first order differential 
equations rather than second order ones. We will show quite generally that all solutions 
of these first order equations preserve (at least) half of the supersymmetries and hence 
are BPS. 

We will first discuss in a general setting when gravity coupled to scalars admits 
solutions of first order equations. Then we will see under which curcumstances this is 
true for the bosonic sector of supergravity, and that in these cases the solutions are BPS. 
Finally, we specialise to the case of present interest and derive the set of coupled first 
order equations of the effective supergravity describing the thermal phases of N = 4 
superstrings. 



3.1 First order equations for gravity coupled to scalars 

Consider the following Z)-dimensional Lagrangian density representing gravity coupled 
to iV scalars 

— C = -R - -Gui^g^d^d^ - V(<p) , (3.1) 

where summation over the repeated indices I, J — 1,2, ... ,N and m,n = 1,2, ... ,D is 
implied. The equations of motion that follow from varying the metric and the scalar 
fields are: 

11 1 

5g : jRmn - -G IJ d m ip I d n ip J = — - ^ g mn V , (3.2) 

5<p : DV + r I JK g mn d m <p J d n <p K = G IJ djV , (3.3) 

where Tj K is the Christoffel symbol formed using the scalar-field space metric Gjj. Let 
us assume that the potential V(ip) is given in terms of a prepotential W(ip) as: 

V = $L {G IJ d I WdjW - 2§^^ 2 ) , (3.4) 

where f3 is a dimensionful constant. 

Let us make the ansatz for a metric of the form 

ds 2 = dr 2 + e 2A{r) r]^dx^dx u , n, v = 1, 2, . . . , D - 1 , (3.5) 
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that preserves the flat space-time symmetries in D — 1 dimensions. We also assume that 
all scalars depend only on the variable r, i.e., if 1 = (p x {r). Then, it can be shown by a 
straightforward calculation that if the following first order equations hold 

the second order equations are also satisfied.^ In proving this statement we use the fact 
that the non-vanishing components of the Riemann tensor for the metric fl3.5p are 



R,u = -v,ue 2A (A" + (D-l)(A') 2 
Rrr = — (D — 1) (A" + (A') 2 ^j , (3.7) 

where the prime denotes the derivative with respect to r. Note that the freedom in the 
choice of sign in (|3.6|) is due to the fact that (|3.4|) does not determine the sign of W. 



Equivalently it expresses the freedom to change r — > — r. 

Throughout the remainder of this paper we restrict to D = 4 and choose (3 = v2 
(consistent with k = \/2). 



3.2 First order equations for the bosonic sector of supergravity 



Consider now the bosonic sector of a generic four- dimensional N = 1 supergravity theory 
as given by eqs. (|2.11 ), ( |2.12 ) and ( p2 . 1 3| ) . Assume that all Im($ / ) vanish and that the 
Kahler potential is such that 



dK 



d 



real directions 



aRe$ 7 



real directions 



(3.8) 



This is indeed the case for our thermal supergravity, but for the time being we do not 
want to consider any specific model. The definition of the covariant derivative fl2.13|) of 
the superpotential W and the assumption (|3.8|) imply that if we define 



W 



,k/2 W 



real directions 



we will have 



_d_ 



W = e K/2 W. I . 



(3.9) 



(3.10) 



As a consequence, the scalar potential (|2.12|) of the supergravity can be rewritten using 
only W as 

1 



V 



_ K adWdw ^ 

4 I d&dQJ , 



(3.11) 



Recently, such first order equations with flat metric in the scalar-field space appeared in the context 
of supersymmetric solutions in five-dimensional gauged supergravity (for a first example, see pq] ). For 
cases where the metric in the scalar-field space is non-trivial, five-dimensional examples have been given 
in |2SH. 
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We see that this is precisely of the form fl3.4j) with D = 4 and (3 = y/2. 



We conclude that if we find solutions to the domain wall ansatz ( |3.5| ) satisfying the 
first order equations 

*i = ±*K"% , £ = , (3.12, 

dr y/2 d$ J dr y/2 v ' 

and take the fermions to vanish, then all equations of motion of the supergravity will be 
satisfied. We now proceed to show that these solutions of ( |3.12| ) precisely preserve half 
or all of the four- dimensional supersymmetries. 



3.3 Supersymmetry of the solution 

It is very simple to obtain the number of supersymmetries left unbroken by a solution of 
the first order equations. The only additional assumptions we need is that the Kahler 
potential satisfies ( |3.8|) , Im($ / ) vanishes, and that the fields only depend on a single 
coordinate, say r, with g rr = 1. 

Consider then the supersymmetry transformation of the (left-handed) fermionic com- 
ponent Xl °f a chiral multiplet with scalar <3> 7 : 

S X [ = -Ly {d^tR - \t L e K/2 K I1 W ]J + fermionic terms . (3.13) 
v 2 2 

Since, by assumption, the solution depends only on a single coordinate r with metric 
g rr = 1, we have 7 At <9 At $ 7 = 7 r <9 r <l )/ where Y satisfies {Y) 2 — 1- Then, using the first 
order equation ( |3.12| ) (and ( |3.10| )) we have 

7 ^„$ / = ±^ 1 r K I1 ^ T = ±-^Ye K/2 K l7 W.j . (3.14) 

' " y/2 d$ J y/2 

It follows that for vanishing fermions 

Sx[ = ~Pl{1 T 7 r )e e^K^W.j , (3.15) 

where Pl is the left-handed chirality projector (of course there is a similar equation for 
Xr}- Clearly, since (7 r ) 2 = 1, the solution will preserve one half of the supersymmetries 
if W-j 7^ 0, and all supersymmetries if W-j = 0, as it should. 

In order to compute the dependence of the Killing spinor e on the space-time co- 
ordinates, one in principle has to solve the gravitino equation. Equivalently, for static 
configurations (as in our domain wall ansatz ( 3.17|) ) one can use an argument based on the 



supersymmetry algebra and deduce that e = (?oo e , where e is a constant spinor subject 
to model-dependent projection conditions, which reduce the number of its independent 
components [00]. Applying this to our case we find 

c(r) = e A ^ 2 e , (lTf)e = 0. (3.16) 
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This result should be interpreted with some care when considering strings at finite tem- 
perature. Formulating a theory at finite temperature requires a rotation to Euclidean 
compact time and, with a trivial background, supersymmetry is completely broken. As 
already discussed in section 2.3, this is the case of the low-temperature phase. Our 
domain wall solutions are, however, backgrounds depending on a spacelike variable r 
with a manifest three-dimensional residual "space-time" symmetry, with signature (3, 0) 
at finite temperature. The supercharges preserved by these backgrounds are then su- 
persymmetries inherited from the original five- dimensional supersymmetry algebra, with 
Minkowski signature (4,1). 

We conclude that the solutions of the first order equations fl3.12p are 1/2 BPS so- 
lutions, except of course if = 0, in which case the scalar fields are constant and 
supersymmetry is unbroken. 



3.4 The first order equations for the thermal supergravity 

Next, we specialise the discussion to the supergravity theory that describes the thermal 
phases of the N = 4 superstrings, by applying the restrictions and identifications ex- 
plained in section 2.2. The Kahler metric, the potential V and prepotential W are given 
in eqs. (|2.21|) , (|2.22f) and (|2.27|) . There are six real scalar fields and the Kahler metric is 



diagonal Kij = K^Sij. For convenience we write the relevant equations again, for D = 4. 
The metric is assumed to be conformally flat and written in the form (see, for instance, 



0) 



ds 2 = dr 2 + e 2A{r) ri^dx^dx v , v = 1, 2, 3 (3.17) 

and the BPS domain wall configurations are solutions of the following system of first 
order differential equations: 

*Pi 1 K -iW . dA 1 w rntt 

taking into account the inverse of the diagonal metric in field space for each one of the 
components </v Note that we have used the freedom r — > — r to make a definite choice 
of signs in these equations. 

The six scalar fields are naturally separated into two groups: {0,; 2 = 1,2,3} corre- 
sponding to s, t, u and the three winding fields {zf, i = 1,2,3}. For the first set the 
metric in field space is flat, i.e., K^((p) = 1, whereas for the second set all components 
are equal but non-trivial, namely Ku)(z) = 4/(1 — z\ — z\ — z 2 ) 2 , see eq. Q2.21| ). The pre- 
potential W was given in ( |2.27| ). As a result, the domain walls of the theory correspond 
to solutions of the non-linear system 



= -^ + ^ + 2e^smh(f )+ H 2 + e~^(e' l) -H 2 2 +e~^H 2 ) , (3.19) 
V2^ = --e^ + ^ + 2e* cosh0 + H 2 - (e^H 2 2 - Hj) , (3.20) 
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V2 



dr 
d in ^ 

dr 

din # 2 





Hi) , 


(3.21) 






(3.22) 


-Hi + e 


-"-Hi) 


, (3.23) 


-Hi + e 


~"-Hl) 


, (3-24) 


= Zi/{1 


-4- 


Z 2 ~ Z l) 



dr 2 1 

^ln#3 = _I e -^- + 4e ^ sinh Hl-2e-* 
dr 2 

where we found it more convenient to work with the fields Hi 
instead of z$ and we defined cj)± = fa ± fo. As soon as a solution has been found, the 
conformal factor of the metric can be obtained by a simple integration of the resulting 
expression for the prepotential W, as it has been prescribed above. 

Although these equations are first order, their explicit solution is still a difficult task 
because they are highly non-linear in the general case. We note that there are examples 
in the scalar sector of maximal gauged supergravity in 4, 5 and 7 space-time dimensions 
where the resulting system of equations can be linearized by introducing appropriate 
combinations of fields plus an algebraic constraint. This led to a systematic classification 
of the possible domain wall solutions in terms of Riemann surfaces with varying genus 
(for details, see the method presented in |32], |3~3l). However, things do not always work 



that way; as it turns out the above system of equations that describe strings at finite 
temperature can only be partially studied by similar methods. 

It can be easily seen from the last three equations that the fields Hi can be solely 
expressed in terms of the fields (pi, but one should avoid substituting their integral ex- 
pressions into the first three equations because it would lead to a complicated system 
of integro-differential equations among the fields (pi alone. In either case, the system of 
equations at hand is rather complicated and difficult to solve in general. Instead, we 
will focus attention on subsectors obtained by consistent truncations of the field content. 
Each one of these consistent truncations results into a system of three first order equa- 
tions that correspond to the various type-II and heterotic theories. Luckily, it will turn 
out that the general solutions for all type-II theories can be found explicitly. However, in 
the heterotic case a general solution cannot be given in closed form. Hence in this case, 
apart from extracting the behaviour of the fields in the strong and weak coupling regions 
as well as around certain critical points, one has unavoidably to rely on numerical stud- 
ies. Combining these pieces of information, we can nevertheless give a rather complete 
picture of how the different solutions behave. 
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4 Type II string theories 

The simplest truncations of the domain wall equations lead to different type-II sectors. 
We will first examine the type-IIA and IIB theories, which can be treated simultaneously 
using the appropriate field identifications, and then examine a hybrid type-II sector, 
which describes type-II strings at the self-dual radius, and turns out to be exactly solvable, 
too. 



4.1 Type-IIA and IIB sector 

According to the field identifications made in the effective supergravity, the type-IIA and 
type-IIB sectors of the theory are obtained by setting 

Zl = Z2 = , for type IIB (4.1) 
Zx = z 3 = , for type IIA (4.2) 

in which case Hi — H< 2 = and Hi = H 3 = respectively. It follows from ( |3.19| ) and 
( |3.20| ) or (|3.21| ) that <pi equals 2 or <p 3 , respectively, up to an irrelevant additive constant 
which we will ignore. Then, it is convenient to set 

\ = 4>2 = ; 03 = x ; z 3 = z for IIB , (4.3) 



or 



V2 



H = 03 = • 2 = x ; z 2 = z for IIA (4.4) 



and introduce in either case the field u for the winding field 

* = tanh(!) , (4.5) 

for which we have 2H = sinhu;. The appropriate string coupling is gn ~ e x and the 
temperature field is T ~ e^^, when it is normalized in type-II string units that remove 
the x-dependence. We treat both cases together because the (pre)potential assumes the 
same form for IIA and IIB, and the same is true for the truncated system of differential 
equations for the domain walls. Hence, no distinction will be made in the following 
between type-IIA or type-IIB. Note also that the kinetic terms for the fields \i an d oj 
all assume their canonical form. 



Explicit calculation shows that in terms of the new variables the prepotential ( |2.27 ) 
becomes 

W n = -e x (e^* + -e-^sinh 2 ^) , (4.6) 
whereas the corresponding potential ( |2.22| ) is 

Vn = -U 2 * sinh 2 oo ( e ~ 2 ^ cosh 2 u - 4) . (4.7) 
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It is clear from this potential Vu that u becomes tachyonic if e~ 2 ^^ < 4 at to = 0, in 
accordance with the discussion of section 2.3. Using the truncated prepotential Wu, the 
type-II domain walls obey the simpler system of first order equations 

^ = \=ex(e^ + L-^smh 2 u 

dr 2y/2 V 2 

d 4 = --e*(e^*- V^siiih 2 u;V (4. 
dr 2 V 2 J ' v 

^ = L e x-^ sinh2 a;. 



These equations can as well be obtained from the full set of first order equations ( |3.19|) 



( p.24 ) by imposing the above field identifications for the type-IIA or IIB strings. The 



equation for the conformal factor of the metric ( 3.17|) is easily integrated and gives A 



— x (up to a constant that can be absorbed into a redefinition of the x^'s). The system 
of equations ( |4.8| ) is invariant under u — > —uo. Therefore we will consider the case u > 
without any loss of generality. 

We note for completeness that in the type-II sector one may fully absorb the \- 
dependence that appears on the right hand side of these differential equations by simply 
changing the coordinate variable r to ( defined as d( = e x dr. Put differently, x plays 
the role of a Liouville field in that the "gravitational dressing" of an auxiliary massive 
model made out of the fields <fr and u 

C((f),u) = ~^(d(j)) 2 -^(duj) 2 -^smh 2 uj(e' 2 ^cosh 2 uj-A) , (4.9) 

yields the scalar field sector of the truncated type-II theory with all three fields 0, u and 
X having kinetic terms in canonical form and a dressed potential that assumes the form 
Vu above. f\ 

It is an interesting finding that the domain wall equations can be completely integrated 
for the type-II sector, and one has for the first time a one-parameter family of explicit 
solutions with non-trivial winding uj. It is convenient for this purpose to present the 
solution by treating the field uo as an independent variable, which is legitimate since the 
third equation in ( |4.8| ) implies that u is a monotonous function of r. Then, the differential 
equation for <p can be easily integrated and also the equation for x- The final result is a 
family of BPS solutions, as shown in section 3.3, with 

e -2V24> = 2 cosh 2 u (In coth 2 u + c) - 2 , 

e 2x = cosh^e^, (4.10) 

parametrized by an arbitrary integration constant c. As we will see, the physical inter- 
pretation of the solutions depends crucially on whether c is positive, negative or zero. 

5 This interpretation is motivated by two-dimensional gravity coupled to massive cr-models. However, 
unlike other examples, where non-abelian Toda theories result in this fashion p4[ , the present model 
does not have a special meaning in integrable systems. 
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The behaviour of the functions e 2 ^ 2 ^ and e 2x for the three choices of c is sketched in 
Figure 1. We also note that we have omitted a multiplicative integration constant on 



the right hand side of the second equation in ( 4.10 ), since it can always be absorbed into 



trivial field redefinitions. We see that, independently of c, the behaviour of T 



~ e 



as u -> is e -2 ^ 
all these solutions. 



21nu;- 



oo. Hence, the winding mode cannot be tachyonic for 




Figure 1: Qualitative behaviour of e 2 ^^ (curves e) and e 2x (curves f) as functions of 
for type-IIA and type- 1 IB. 



Performing the integration of the conformal factor, we find that the metric takes the 
form 

ds 2 = o i — du) 2 H ^ — Vuvdx^dx" . (4-H) 

sinh cucosh lo cosh lo 

As for the relation between the variables r and u in ( p. 17] ) and ( [4.11 ) this turns out to 

be given by the relation of differentials 

dr = 2y ^ eH 2 du . (4.12) 

sinh u cosh u 

Clearly, this integration cannot be performed in closed form and so we lack the explicit 
dependence of u(r). However, the dependence can easily be spelled out in some limiting 
cases to which we now turn. 
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Universal behaviour for small windings: First consider the limit of vanishing winding 
field to. Then the asymptotic behaviour of the fields is 

/ o \2/3 

w ~ e -(3r/8)^ ) e -2 X „ e -V2^ ^ M r ; as r ^ oo . (4.13) 

Hence, in the limit of vanishing winding field uj, the fields <fi and \ approach — oo and 
the value of the constant c plays no role. [] In addition, the metric becomes 

( 3 V /3 

ds 2 ~ dr 2 + —j= r r\ ilv dx p 'dx v , as r^oo. (4-14) 



The behaviour of the solution for larger values of the winding uj does depend on the 
parameter c. We distinguish the following three cases: 

c > : In this case the reality conditions for the fields x an d 4> allow for < uj < oo. 
Then, the string coupling e 2x is a monotonically increasing function of uj in its full range 
of values. However, the temperature field squared e 2 ^ 2 ^ first increases until it reaches a 
maximum at uj = ujq and then decreases to zero as uj — > oo. In the limits of small and 
large c the constant uj can be found analytically: 

f —I In ( 2.) for c -> 0+ 
^o= / W ' (4.15) 

I for c —y oo . 

For intermediate values of c, c<Jo ranges between the above two extremes. The maximum 
value of e 2 ^ 2 ^ is given by 



e 2v / 20 



= -sinh 2 ^ . (4.16) 

OJ=Ll>0 2 



The asymptotic behaviour of the various fields is found by noting that due to ( 4. 12] ) we 
have the relation 

32 • 2 3/4 

r ~ e~ 7ul/2 , as r^0 + and w^oo, (4.17) 

up to an additive constant that has been fixed in an obvious way. Then, 

e-^ ~ (^jP) ' , z~ 2x - (^ cV ) 2A ' as r -> 0+ . (4.18) 
Similarly, the asymptotic expression for the metric is 

( 7 \ 2A 

ds 2 ~ dr 2 + l-j=c 2 rj rj^dx^dx" , as r -> + . (4.19) 

c = : As before, the reality conditions for the fields x an d allow for < uj < oo. 
However, both fields e 2x and e 2 ^ 2 ^ are now monotonically increasing functions of uj in 

6 The exponential vanishing of the winding field uj is related to the fact that setting uj — from the 
very beginning is consistent with the system of eqs. (4.8). In fact, then the solution ( |4 .13 ) for \ an d <t> 



and the metric in ( |4. 14J) below become exact (see also (|5.7| ) below) 
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the entire range of values. The asymptotic behaviour of the various fields is found by 
noting that due to ( |4.12[ ) we have the relation 

32 



Then, 



r ~ — e 5uj/2 , as r^0 + and w^oo. (4.20) 
5 

e -V2^ _ i( 5r )2/5 ? e -2 X _ I( 5r )6/5 ? as r ^ 0+ . (4.21) 

2 8 
Similarly, the asymptotic expression for the metric is 

ds 2 ~ dr 2 + l(5r) 6/5 ri^dx^dx" , as r -> + . (4.22) 
8 

c < : In this case the reality conditions for the fields % an d do not allow to exceed a 
maximum value for the winding field, i.e., < uj < u> max . In this range, e 2x and e 2 ^ 2 ^ 
are monotonically increasing functions of uj. Similarly, in the limit of small or large |c| 
the constant cj max can be found analytically, as before, 

f r < ~°~' (4.23) 
lea for c — > — oo , 

whereas for intermediate values of c, tu max ranges between the two above extremes. The 
asymptotic behaviour of the various fields can be found by noting that due to ( |4.12| ) we 
have the relation 



o 3/4 7/ 4 V~max 
OSm Cm 



uf'\ as r^0+ and u -+ < ax , (4.24) 



by introducing the notation 

s m = sinhu; max , c m = coshu; max . (4.25) 



Then, 



-^(^r) , ^J' r] , as _0 + . (4,6, 




Similarly, the asymptotic expression for the metric is 

/ 3 \ 2/3 

ds 2 ~ dr 2 + — -= — r rj^dx^dx" , as r -> + . (4.27) 



In all cases above, we note that the main qualitative difference is in the shape of 
the function e 2 ^^ representing the temperature field square: for c < it starts from 
zero and evolves monotonically to +oo following the range of r from oo to 0, whereas 
for c > we see a dramatic change in that the function starts and ends at zero, thus 
developing a maximum along the way. The end point at r = corresponds to a curvature 
singularity. As we will see in section 6, there are physical ways to differentiate among the 
allowed values of c by requiring consistent propagation of a quantum test particle on the 
corresponding supergravity backgrounds which have a curvature singularity. Of course, 
string theory (and not supergravity alone) holds the ultimate answer for the physical 
relevance of our solutions in thermodynamics. 
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4.2 A hybrid type-II sector 



There is another truncation of the domain wall equations with H\ = 0, as in type-II, 
which is exactly solvable. We set for this purpose H 2 = ±#3 and observe that the full 
system of equations is consistent provided that 02 = 03- Actually, this sector can be 
viewed as a hybrid of type-IIB and IIA in that we impose the "diagonal" constraint 
H 2 = ±i^3 instead of choosing the axes H 2 = or H3 = respectively in if-field space; 
as such, it is a hybrid truncation. This truncation amounts to choosing the self-dual 
radius so that there is no dinstinction between the type-IIA and type-IIB theories. 

We proceed further by setting 

0i = X , 02 = 03 = ^0 , H 2 = ±H 3 = ^|Sinh^ (4.28) 



and so z 2 = ±z 3 = z can be chosen as 



^tanh| , (4.29) 



in order for the kinetic terms of the fields Xi an d oj to assume their canonical form. As 
before, the temperature field is T ~ e^^. Then, the prepotential (|2.27 ) and the potential 
( PI become 

W hyh = - (V +v/ ^ + ^e-*sinh 2 ^ (4.30) 

and 

Vhyb = — sinh 2 uj (e- 2x {\ + cosh 2 u) - 8e^) , (4.31) 
32 

respectively. We see that u will become tachyonic, if at uj = we have er^^~ 2x < 4. 
The truncated system of equations is 



dr 2\/2 V 2 J 



2V2 

—— = =e~ x sinh 2uj . 

dr 4V2 



The equation for the conformal factor of the metric ([3.17 ) can be easily integrated; it 



gives A = ~x — ln(cosha;), up to a constant that can be absorbed into a redefinition 
of the x M 's. Notice that the result is not the same as in the genuine type-II case that 
we examined before. The system of equations ( |4.32|) is invariant under uj — > — 00, and as 



before we only need to consider the case of u > 0. Note also that, unlike the genuine 
type-II case, the field x cannot be viewed as a Liouville field in this hybrid sector that 
provides the "gravitational dressing" of a simpler massive model for the scalar fields 
and uj. 
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Luckily, the general solution can be found in parametric form by employing uo as an 
independent variable, as before. We find the following family of BPS solutions 



= 6 + 2a [In 



coth^ 



1 



cosher 



e 2 X = 



. , (4.33) 
2coshu> 

where a and b are integration constants with a being positive for reality reasons. For 
convenience, in the rest of this subsection, we set the positive constant a = 2, since its 
precise value can be adjusted by field redefinitions and hence it has no physical relevance. 
The other integration constant b cannot be fixed mathematically, but we will see later that 
physical considerations in the context of supergravity may impose some restrictions on 
it. Of course, as before, it is string theory that holds the answer for having an acceptable 
geometrical background for domain walls. The behaviour of the functions e 2 ^* and e 2x 
for the three choices of b is sketched in Figure 2. Again, as u> — > 0, the temperature field 
behaves as T' 2 ~ e - 2 ^<t> ~ 21ncu~ 2 — > oo and e~ 2x ~ — > oo so that cu cannot be 

tachyonic. 




Figure 2: Qualitative behaviour of e 2 ^ 2 ^ (curves e) and e 2x (curves f) as functions of uo 
for the hybrid type-II. 



After computing the conformal factor we may write the metric as 



ds 2 



sinh cjcosh u 



2 e -v / 20 

dui 2 H — rjnudx^dx" . 

a cosh uj 



(4.34) 
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The connection between the variables r and to in ([3.17 ) and ( f4.34| ) is given by the relation 
of differentials 

2Jae <t>/ ^ 2 

dr = — duo . (4.35) 

sinh uj cosh ' uo 

As before, the integration cannot be performed in closed form, apart from a few limiting 
cases, as we will see next. 

Universal behaviour for small windings: Let us consider first the case of vanishing wind- 
ing field uj. It turns out that the asymptotic behaviour of the fields is 

/ o \2/3 

uj ~ 2e~ {3r ^ /3 , e~ 2x ~ ~ I -= r) , as r -> oo . (4.36) 



Hence, in the limit of vanishing winding field uj, the fields <fi and \ approach — oo and 
the value of the constant b plays no role.[] The metric becomes 

/ 3 \ 2/3 

ds 2 ~ dr 2 + —= r r/^dx^dx 1 " , as r^oo. (4.37) 



Note that this asymptotic behaviour is almost identical to the one found for the pure 
type-II case. A similar behaviour will also exist in the heterotic case to be discussed 
in the next section. There we will see that it is related to an exact solution with zero 
winding present in all cases. 

The behaviour of the solution for larger values of the winding uj, however, depends 
on the parameter b. As in the genuine type-II solution ( |4.10|) , we also distinguish here 
the following three cases: 

b > : In this case the reality conditions for the fields x an d 4> allow for < uj < oo. 
Then, the temperature field square e 2 ^^ is a monotonically increasing function of uj until 
it reaches the asymptotic value 1/6. However, the string coupling e 2x first increases until 
it reaches a maximum at uj = ujq, and then decreases to zero as uj —>■ oo. In the limits of 
small and large b, the constant uj q can be found analytically 



-| In (f|) for 6^0- 
for b — » oo 



^o=<| 1 U6y r , ' ' (4.38) 



whereas for intermediate values of b, ujq ranges between the above two extremes. The 
maximum value of e 2x is given by 



e 2 * 



sinh ujq 

(4.39) 



\J2 cosh uj 



The asymptotic behaviour of the various fields is extracted by noting that, due to ( |4.35| ), 
we have the relation 



32 

— Tje- 5 ^ 2 , as r^0 + and w^co. (4.40) 
5o i/4 



7 A similar comment, as in footnote 6, applies here as well. 
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Then, we find that 

e -2V20 „ b + 1 ( 56 l/4 r) 6/fi ? e 2 X „ 1 ^^) 2/5 ^ as r ^ + (4 41) 

Similarly, the asymptotic expression for the metric is 

ds 2 ~ dr 2 + ^(56 3/2 r) 2/5 ^rf^cix y , as r -»• 0+ . (4.42) 

6 = : As before, the reality conditions for the fields \ and <fi allow for < to < oo. 
However, both fields e 2x and e 2 ^ 2 ^ are now monotonically increasing functions of u in 
the entire range of values. The asymptotic behaviour of the various fields is found using 

' ' e- 7uj/4 , as r^0 + and w-^oo, (4.43) 



7 V2 

which is due to ( 4.35| ). Then, we find that 



e -2v^ „ 7_ V7 ^ ^ ^ 2x ^ ^|8^ l ' 7 r2/7 ^ ^ r ^ Q+ (4 M) 

Similarly, the asymptotic expression for the metric is 

7 /3-7 3 \ 1/7 

rfs 2 ~ dr 2 + — f — ^— J r w/7 r]^dx^dx u , as r -> 0+ . (4.45) 

b < : In this case the reality conditions for the fields x and do not allow to exceed a 
maximum value for the winding field, i.e., < u < u; max . In this range, e 2x and e 2 ^ 2 ^ 
are monotonically increasing functions of uj. In the limit of small and large |6| we find 
for c<j max the analytic expressions 

f -Ilnf^) for 0" , 

w max - I * A } 32 * 4.46 
\ 2e fe / 4 - 1 for 6 -> -oo . V ; 

For intermediate values of 6, w ma x ranges between the two above extremes. The asymp- 
totic behaviour of the various fields is extracted by noting that, due to (|4.35|) , we have 
the relation 



3s 3/4 c ' 

where we have used the same notation as in fl4.25|) . Then 



3/4 , as r -> + and u -> w max , (4.47) 



e -^~ f— ^=r") /3 , e- 2 *~(%rV /3 , as r ^ 0+ . (4.48) 

^v 7 ^ J {2V2 J K ' 

Similarly, the asymptotic expression for the metric is 

/ 3 \ 2/3 

ds 2 ~ dr 2 + f -^-j= - rj r] fl „dx fl dx u , as r — > + . (4.49) 
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Note that in the hybrid type-II sector, as for the genuine type-II solution, the end 
point at r = represents a curvature singularity. We will see in section 6 that the three 
different cases, which correspond to the values of b > 0, or < 0, can be distinguished 
using certain physical criteria. 

Finally, note that for all type-II sectors (IIA, IIB and the hybrid type-II), u decreases 
monotonically from its value at r = (which may be infinite) to u = at r = oo. 
In particular, this implies that uo{r) cannot be a periodic function of r, and hence r 
necessarily is a non-compact coordinate. 



5 Heterotic sector 

The heterotic limit is obtained by setting H2 = H3 = 0, while keeping H\ = H free to 
vary. It then follows from (|3.20|) and ( |3.21| ) that 02 = 03 up to an irrelevant additive 



constant that we will ignore. Trading the z corresponding to H for a winding field u and 
introducing appropriate variables 

0i = X, 02 = 03 = ^, z = tanh^, (5.1) 

for which 2H = sinhcu, we obtain the expression for the truncated prepotential 

W het = ^e x (e^ - sinh(V20)sinh 2 o;) , (5.2) 

and the potential 

V hct = -e 2x sinh 2 u (cosh(2 v / 20) cosh 2 u - 3) . (5.3) 
8 ^ ' 

We note for completeness that e x is the string coupling and the temperature field is 
again (normalized in heterotic string units that remove the x-dependence) and the 
kinetic terms of the fields Xi an d oj assume their canonical form with our choice of 
normalizations. We see that to can be tachyonic only if cosh(2y / 20) < 3 at to = 0. 

We remark that although the potential Vhet is invariant under uo — > — u and — > —0, 
the prepotential Whet is invariant under the first transformation only. The second trans- 
formation corresponds to T-duality in the temperature radius. Since the prepotential 
breaks this temperature T-duality one expects the same to be true for non-trivial solu- 
tions of our equations. 

The truncated system of differential equations is 

^ = 1 - r e x (e^ - sinh(v^0)sinhV) , 

(XT 2 v 2 

l e x (eV24> _ cos h(v / 20)sinh 2 cj) , (5.4) 
re x sinh(\/20) sinh 2a; . 



dr 2 
du 1 



dr 2V2 
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The equation for the conformal factor of the metric (|3.17|) is easily integrated and gives 
A = ~x, as in the case of the genuine type-IIA or type-IIB theories. Also, similarly to 
this case, we may drop the x-dependence from the right hand side of the resulting first 
order system by simply changing variables to ( defined by d( = e x dr. Again, x may 
be viewed as a Liouville field whose coupling provides the "gravitational dressing" of an 
auxiliary massive model of the fields and uj with Lagrangian density 

£'((/), u) = --(90) 2 --(^) 2 --sinh 2 w(cosh(2 v / 20)cosh 2 cj-3) . (5.5) 

This prescription produces naturally the canonical kinetic term for the field \ and dresses 
the potential into Vh et . 

It turns out that it is not possible to solve the domain wall equations in closed form in 
the heterotic sector. Nevertheless, we can still determine the leading small to asymptotics, 
similarly to the analogous universal behaviour for the type-II theories. This corresponds 
to large r. As it turns out, to is again exponentially suppressed as r — > oo, while e x and 
gv 7 ^ h ave a power-law behaviour. Hence, we can safely drop the sinh 2 uj part in the first 
two equations in (|5.4j ) and replace sinh2co> by 2uj in the third equation. The resulting 
system is easily integrated and yields in this limit 

e -V5* „ 2CV 2/3 ^ e -2 X „ 9 r2/3 ^ ^ ^ g _ c2r 4/3 ^ as r ^ oo , (5. 6 ) 

where C and ujq are two positive constants of integration. Note that cosh(2\/20) — 
2C 2 r 4//3 ~ 21na;~ 2 — ► oo as uj — > and uj cannot be tachyonic. As one decreases the 
value of the constant ujq, the range of approximate validity of this asymptotic solution 
extends to larger intervals of r. Eventually, if ujq = one obtains a particular exact 
solution^] 

e-^ = 2Cr 2 / 3 , e- 2 * = ^r 2 / 3 , w = 0, (5.7) 

which is valid for all values of r. However, since we lack the most general solution, even in 
some parametric form <f)(uj), we will only provide results about the dominant behaviour 
of the fields in the weak coupling (i.e., x ~ > — oo) and the strong coupling (i.e., x ~^ +oo) 
regions. In fact, we will succeed to identify these regions and perform an in depth analysis 
in their vicinity using standard techniques from dynamical systems. The remaining part 
of the two-dimensional space (0, uj) can only be studied numerically. 

Before delving into details, we study a bit more the exact solution for uj = 0. This 
solution represents a straight orbit in the parameter space (0, uj) with uj = everywhere, 
extending from = — oo (identified as a weak coupling point) to = +oo (identified as 
a strong coupling point), as r ranges from +oo to 0. The Liouville coordinate (, which 
trades r as d( = e x dr, is given for this special solution by 

C = 2 v / 2Cr 2/3 , (5.8) 

8 We note that this is an exact solution for the type-II sector as well. Indeed, when uj = 0, the three 
different truncated systems of equations ( |4.S| ), ( I.32| ) and (5.4) become identical. Note that for the 



corresponding type-II solution we had chosen the integration constant C = 3 2 / 3 /4. 
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where we absorbed an integration constant into £. 

Introducing small fluctuations 50 and 5uj around this special solution, to test its 
stability, we find to first order that they satisfy the linearized system of equations 

dC = 

^5u = [--^)5u, (5.9) 

which are written for convenience using (. Therefore, we have the following solution of 
the linearized problem 

<ty=4> 6uJ = B \f(e~ <2/8 , (5-10) 

where A, B are some other constants. These fluctuations are small provided that ( — > oo, 
otherwise the linearized approximation is not valid. We conclude this analysis by noting 
that the exact solution we found is stable against perturbations in the parameter space 
(<p,co) provided that we are in the weak coupling region where — * — oo. In fact, as we 
will see next, the point (—00, 0) is an asymptotic critical point in the two-dimensional 
parameter space, and so we can talk about asymptotic stability in its vicinity. 

The critical points: It is more convenient to work with the coordinate ( instead of r in 
order to drop the x-dependence from the right hand side of the equations. Then, as in 
the theory of dynamical systems, we have to study the solutions of the system 

^ = P(0,w), ^ = W,w), (5.11) 

either as functions of (, i.e., as <p(() and or in a parametric form <p{uj). For the 

moment we have not prescribed a set of physical boundary conditions that select certain 
orbits among the infinite many that arise in the two-dimensional parameter space (0, u), 
as this will be done in the next section. Here, we only focus on general properties of the 
solutions around the critical points of the system. 

Recall that the critical points are defined as the common zeros of the functions P and 
Q. These occur in the present case either when = and sinhtu = ±1 or when = 
and to — 0. The second case has already been analyzed and lies in the weak coupling 
region of the exact solution that was described above with approaching asymptotically 
—00. Therefore, we focus on the other critical points 



O = O, = ln(V2 ± 1) , (5.12) 

which appear symmetrically on the tu-axis because of the invariance of the equations 
under to — ► —to. Then, the dilaton equation in ( |5.4j ) is easily integrated to give at these 
critical points 
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where we have absorbed integration constants into £ and r in an obvious way. Hence, the 
dilaton field exhibits a logarithmic dependence on r and a linear in (. Recall now that 
the string frame metric in four dimensions is obtained by ( 3.1 7|) after multiplying with 



e x . Hence, the resulting backgound has a flat string metric and a linear dilaton, which is 
the exact solution preserving half of the supersymmetries found earlier in [IB, Note 
that, in this case, it makes sense to pass from the Einstein to the string frame since, 
when and u assume their critical values ( |5.12| ), we are left with only massless fields 
(graviton and dilaton) that couple to a 2-dimensional string world-sheet action. 

Next, if we linearize around these critical points ( |5.12| ) by introducing small fluctua- 
tions 8(f) and Su, we will find the following evolution for them 



d 



f 50 ^ 


1 


( 1 T2\ 




( W \ 


V Su J 


~ V2 


KT2 J 




\ fa ) 



(5.14) 



whereas for the dilaton d(5x)/d( = 0. Computing the two eigenvalues of the correspond- 
ing matrix we have in either case of uJq that 



A- 



2^/2 



-1 ± V17) 



(5.15) 



and since A + > and A_ < 0, we see that the pair of critical points (0,Uq) are both 
saddle points. We may solve the linearized system of equations and obtain 



5(f) = Ae x+<: + Be x - C 
' A 



±v ~ 2 r 



B 
AT' 



>A-C 



(5.16) 



with ± depending on the choice of critical point uJq, whereas the dilaton does not change, 
to linear order, from its value in ( |5.13| ). 

If we demand that the critical points are reached as ( — > +oo, we conclude for the 
integration constants that A = 0. Hence, the variations of <fi and u have a power law 
decay, namely 



VT7 + 1, 
=F ouj 



B 



17+1 



as r 



(5.17) 



On the other hand, if we demand instead, that the critical points are reached as ( — ► — 00, 
we conclude for the integration constants that B = 0. Now eq. ( |5.17|) is replaced by 



17- 1 



717-1 



5u = A 



2V2, 



as r 



(5.18) 



Note the different behaviour of the two solutions around the critical points ( |5.12|) as 
given by ( |5.17p and (|5.18| ). In particular, from ( |5.13|) , they correspond to weak and 
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region 1 : 


0" 


-> +oo 


and 


uj - 




region 2 : 


<t>~ 


-> +oo 


and 


UJ - 


-> ±oo 


region 3 : 


0" 


-> — oo 


and 


UJ - 


-> ±oo 



strong string couplings, respectively. We will see that if trajectories in the entire <fi-uj 
plane are considered, both solutions will appear. 

Strong coupling points: Having extracted the behaviour of the fields around the critical 
points, which can be either at weak or at strong coupling, we will examine next their 
behaviour in the vicinity of the other strong coupling points, where the dilaton field x 
tends to +oo. We have identified three such regions in the parameter space (0, uj), which 
of course have a mirror counterpart under uj — > — uj. They are (cf. figure 3 below): 

(5.19) 
(5.20) 
(5.21) 

To justify this claim and extract the details of the fields in these regions, we shall treat 
each case separately. 

First, it is convenient to examine the heterotic equations in the limit (f> — > +oo, which 
is common to both regions 1 and 2. Using the coordinate ( (instead of r) we obtain in 
this limit the simplified system 

^ = — ^e^sinh^coshu; (5.22) 
«C 2y 2 

and so we see immediately that a/2% — (up to an additive constant). The remaining 
equations can be easily integrated to yield the expressions 



2a^2 



(6-aC) 1/3 (l-(&-aC) 2/3 ), 



coshV, = {b _ ^ c)2/3 , (5.23) 

in terms of some constants a > and b. However, recall at this point that we are 
considering the behaviour of the equations for <fi — > +oo, and so we should expand 
exp(— v^20) around zero. Naturally, we face two possibilities, which will be identified 
with regions 1 and 2. 

Region 1 follows in this context by assuming (b — aQ 2 ^ ~ 1. Defining the parameter 
£i = (b — I) /a we have 

* ^(C " Ci) , a; ~ ±^|(C - Ci) 1/2 5 as C - Ci + (5.24) 

and so uj — > as advertized for region 1. Since V2x = 4> ( U P to a constant shift), we 
conclude that it is indeed a strong coupling region. One easily sees that (C — Ci) ~ 
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(r — rx) 2 / 3 , for some constant r\ which we henceforth absorb into r so that the metric 
behaves like 

ds 2 ~ dr 2 + (const.) r 2/3 r] fll/ dx tJ, dx u , as r — >■ + . (5.25) 

Region 2 follows by expanding the fields around b — a( ~ 0. If we let C2 = b/a, we 
find the following behaviour 

In this vicinity we have ui — * ±00, depending on the branch in the upper or lower a; 
half-plane, and since <fi — > 00 we also have x - > +00, which indeed describes the strong 
coupling region 2 as advertized. Using (£ 2 — C) ~ ( r 2 — r) 6 / 7 , for some positive r 2 , one 
finds that 

ds 2 ~ <ir 2 + (const.) (r 2 - r) 2 ! 1 r} iliV dx ll dx v , as r —> r 2 . (5.27) 

Region 3 arises by first letting <fi — > — 00. Then, the heterotic equations simplify in 
this limit to the following system: 

# = -72^ = -e-^sinh 2 ^ , 

^ = ^e-^sinhcucoshcj . (5.28) 

dQ 2y 2 

Clearly, v^x — — 4> ( U P to a constant) and so this is also a strong coupling region. Since 
we have <p — > —00, the dominant dependence on ( is easy to extract. We find 



1/3 /o./^/s 



e 



2V2\ /d 1 



2c^J (C " C3) ' 6 "H^J (C^' aSC ^ C3+ ' (5 ' 29) 



where c and £3 are constants of integration. In this case we see that — > ±00 as it was 
initially stated. Finally, £ — £3 ~ (r — r 3 ) 6 / 7 for some constant r 3 which we again absorb 
into r. The behaviour of the metric near r = is 

ds 2 ~ dr 2 + (const.) r 2/ ' 'rj^dx^dx" , as r -> + . (5.30) 

The behaviour of different trajectories in the (0, u;) plane is depicted in Figure 3. 
Because of the symmetry uj — > — uj we have restricted to windings uj > without loss of 
generality, in particular solutions that start in the upper-half plane cannot end somewhere 
in the lower half-plane. The qualitative characteristics of the various trajectories can be 
deduced by treating uj (or 0) as an independent variable and consider the differential 
equation for d<p/duj. Before we examine the different trajectories in some detail, note 
by inspection of the differential equation dx/duo and from the strong coupling behaviour 
that we already know, that there can be trajectories along which the string coupling e x 
develops a minimum at a point in the (p-ui plane where 

e 2V2<t> = s^ 2 " ^ $ > o and u > + = ln(v /3 + (5 31) 

sinh a; — 2 
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Region 3 




Region 2 



Region 1 



Figure 3: Qualitative behaviour of different trajectories in the upper-half <p-u plane for 
the heterotic case. 



For <p < there is no such minimum and the string coupling e x keeps increasing with 
increasing uj. Let us now discuss the different types of trajectories as depicted in Figure 
3: 



Trajectories with an extreme uj: These are trajectories where for some uj we have 







and hence duj/dip 
6 — 0. Since 



showing that these trajectories intersect perpendicularly the axis 
d 2 uj sinh 2uj 



0=0 



sinh 2 uj — 1 ' 



(5.32) 



o 



that this uj is a minimum (maximum) in 



we have for uj > uJq = \n(\/2 + 1 
the trajectory. Two such typical trajectories are drawn in Figure 3. The lower one 
connects the weak coupling region (<f),uj) = (— oo, 0) with the strong coupling region 1, 
after reaching a maximum for the winding uJ maiX G (0,uJq). Since ( |5.31|) can never be 
satisfied, it can be easily seen that dx/d<p > throughout the trajectory and hence the 
string coupling e x is progressively increasing with 0. Another such trajectory connects 
the strong coupling regions 2 and 3 going through a minimum for the winding co> m i n G 
(uJq, oo). For this trajectory there is of course a minimum for the string coupling at a 
point described by ( |5.31| ). We also note that for the trajectories with uj < u maiX that 
connect the weak coupling region to region 1, the variable ( and hence also r decrease as 
(f) increases, with r going from +oo to 0. For the trajectories with uj > uj m i n that connect 
regions 3 and 2, ( and r increase as <p increases, with r going from to some r2 > 0. 
Note that the value of r 2 depends on the chosen trajectory. 

Trajectories with an extreme <p: There are trajectories for which d(p/duj = and hence 
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there is an extreme value for 0. We find that the condition for this is 



e 2V24> = smh 2 g ^ | > , for w + <w<w^, ^ ^ 



2 — sinh uj \ < , for 0<w<w 



Since 



du 2 



V2 



sinh u> — 1 



(5.34) 



we have for > that this is a minimum (indicated by m i n ) occuring for u G (wq", UJ\), 
whereas for < it is a maximum (indicated by max ) occuring for uj G (0, Wq"). Two 
such typical trajectories are drawn in Figure 3. The one with < connects the weakly 
coupled region having (0, uj) = (— oo,0) with the strongly coupled region 3. Since < 
there is no mimimum for the string coupling, which increases with u in accordance with 
( 5.31] ). The trajectory with > connects the two strongly coupled regions 1 and 2. 



For such trajectories, besides the depicted minimum value for 0, there is also a minimum 
value for the string coupling according to (|5.31| ). Note that for the trajectory with > 
(0 < 0) the variables ( and r increase (decrease) as uj increases, with r ranging from 
to some T2 > (from oo to 0). Again, the value of r<i depends on the chosen trajectory. 

We also note that region 2 contains trajectories that come either from region 1 or from 
region 3.0 The family of such trajectories is parametrized by the constant a in ( 5.26Q . 



The values of T2 and a are clearly related, but we do not know exactly how. There exists 
a critical value, which is numerically found to be a cr i t — 0.9, such that for a > a cv \t 
(a < a^it) we connect to region 1 (3). Similarly, in region 3 the different trajectories 
are parametrized by the constant c in (|5.29|) and there is also a critical value c CT \ t — 1.2; 
for c > c crit (c < c crit ) these trajectories end up in the weakly coupled region (region 
2). Of course, for the trajectories that connects the regions 2 and 3 the constants a and 
c are related, but how precisely, cannot be answered without knowledge of the explicit 
solution. 

Trajectories ending at critical points: These are the trajectories depicted with the dashed 
lines and occur when the constants a and c assume the critical values we mentioned. 
We emphasize that despite appearances, these trajectories have the critical points as 
their end-points, and they do not simply pass through them. In other words, they 
do not connect the strongly coupled regions 1 and 3 nor the weakly coupled region 
(0,o;) = (—oo,0) with the strongly coupled region 2. The critical point indicated in 
Figure 3, represents a weak coupling point for trajectories that come either from region 
1 or from region 3. In these cases the solution in the vicinity of the critical point is given 
by (|5.171) and the variables ( or r increase to +00 as we approach the critical point. 
In contrast, for trajectories that go either to region 2 or to the weakly coupled region 
with (0,o;) = (—oo,0), this critical point is at strong coupling. Then, the solution in 
its vicinity is given by (|5.18 ) and the variables ( or r increase as we get away from the 



9 In these regions we have the linear behaviour ^/2<f) ~ ui + \xi(^/2a/?>) (for region 2) and \/2<p 
-lo + In (2/c) (for region 3). For region 1 we have instead V2(f> ~ —2 \nui — ln(3/2\/2a) . 
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critical point. We note that the trajectory going to region 2 connects strong coupling to 
strong coupling so that \ must have a minimum along its way. This is consistent with 
the fact that such a minimum occurs if and only if > and u > uj^. This trajectory 
is the only one among the four for which r takes values in a finite interval [0, r 2 ]. 

In section 6, we will discuss certain criteria based on the propagation of a free quan- 
tum particle on the corresponding backgrounds that render some of these trajectories 
unphysical. 

We conclude the discussion of the heterotic equations by mentioning the change of 
variables 

, r i\ sinhu; 

x = tanhu; , y = cosh(V2<p) -5— , (5.35) 

cosh uj 

which cast the differential equations for 4>{ui) into a simpler looking form for y(x) (after 
taking its square) 

x 2 (y') 2 - y 2 = x 2 (x 2 - 1) . (5.36) 

Unfortunately, we are still lacking its general solution in analytic form. Nevertheless, 
the form of the eq. (|5.36|) enables to perform a systematic study of the perturbative 



expansion around the critical points (|5.12 ), which are located at (x, y) = ± (4j, |) in 



the parametrization Q5.35|) . We have the following perturbative expansion, with the upper 



(lower) signs refering to the corresponding signs in ( |5.12j ), 



n=0 

where the coefficient ao of the expansion is given by 

.2 



±0 + E a n [x T ^= j , (5-37) 



2a 2 T a - 2 = . (5.38) 
We also have the following recursive relations for n > 1: 



2 n— 1 

[2(n + 2)a Tl)a n = ±2\plb n ^ - <5 n>2 - - ^2(k + 2)(n - k + 2)a k a. 

1 k=i 



n-l 



TV2j2( k + 2 )( n - k + 1 ) a k a n-k-i (5.39) 

k=0 

n-2 

~ E(( fc + 2 )( n l)a k a n - k - 2 ■ 



k=0 



Note that the freedom to choose the constant ao in (|5.38|) corresponds to the two inde- 



pendent solutions of (|5.36|) that arise from taking the square root. Since the constant a 



determines the sign of y', we have for uj > that ao > for <p > 0, whereas ao < for 
<f) < 0. The first few terms of the expansion are given by 

1 / 1 \ 2 2V2(a ±l) / 1 \ 3 
y = ± - +ao ( XT - )-—L—li XT -\ +.... (5.40) 
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One could develop a similar perturbative expansion around any of the strongly coupled 
regions 1, 2 and 3 or in the vicinity of the weakly coupled region (— oo,0). The corre- 
sponding points, in the x-y plane, around which one should expand are (x, y) = (0, +oo), 
(l,v^a/3), (l,c/2) and (0,0) respectively 

Finally, let us comment on the possibility to have a compact coordinate r as suggested 
by the microcanonical analysis of refs. JOI [16], O. This would imply that all fields must 
come back to their initial value after some period in r. We have seen that there are 
two types of trajectories. On the one hand, we have those for which r takes values in a 
semi-infinite interval. These are the trajectories that connect the weak-coupling region 
{(f), ui) = (—oo,0) to any other region, as well as those connecting regions 3 or 1 to the 
critical point. Clearly, these solutions are not, and cannot be made periodic in r. On 
the other hand, we have the trajectories ending in region 2 for which r takes values in 
the finite intervals [0,r 2 ] with the value of r 2 depending on the chosen trajectory. This 
is crucially different from what we found in the type-II cases. Note that trajectories in 
this second class always go from strong to strong coupling. It may well be that these 
solutions can be continued into periodic functions of r beyond these intervals, tracing the 
corresponding trajectories back and forth. The situation is examplified by the parabola 
y = x 2 in the x — y plane with x a function of r given by x(r) = In ^osr • This clearly is 
a periodic function in r and hence r can be taken as a compact coordinate with values 
on the circle [0,27r]. On the other hand, if we only look at the interval r e [0, 7r], x(r) 
goes from — oo to +oo and we trace the parabola just once. This is somewhat similar 
to our trajectories going into region 2. Without an explicit solution at our disposal, it 
is probably impossible to decide whether or not we can continue these solutions into 
periodic ones or not. 



6 Physical boundary conditions 

Next, we examine the boundary conditions which are physically relevant for the domain 
wall backgrounds. Note that until now the integration constants of the corresponding 
first order equations were left undetermined, and so any solution could be useful math- 
ematically, as there is no preference among the different orbits in the parameter space, 
say (<f), u) . We now apply some general criteria for restricting the physical range of the 
parameters within the context of supergravity using the propagation of a quantum test 
particle on the corresponding backgrounds. Then, we discuss the criteria that string 
theory may impose, although the situation is less clear in that case. 

Curvature singularities lead to singular classical dynamics for test particles. If they 
persist at the quantum level, the theory is considered as unphysical. Unfortunately, 
we do not know enough about string theory in such backgrounds in order to answer this 
question definitely. Instead, we can ask the simpler question whether the propagation of a 
quantum test particle is well defined in the presence of a curvature singularity. A singular 
classical propagation, as indicated by an incomplete geodesic motion, does not necessarily 
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lead to a singular wave propagation since the space-time can produce an effective barrier 
that shields the classical singularity. In general terms, a singularity will be physical if 
the evolution of any state is uniquely defined for all times. The background is then 
called wave-regular. This criterion is identical to finding a unique self-adjoint extension 
of the wave-operator |35], [56], [37] . This can be quickly investigated by looking at the two 



solutions of the wave equation locally near the singularity. Note that the wave equation 
arises, for instance, when one considers linearized graviton fluctuations along the three- 
dimensional brane embedded into the four- dimensional domain wall ansatz ( |3.17| ), i.e., 



Vnu ~~ * V^f + ^m^-Q m sucn cases 5 m addition to having a unique evolution of initial data, 
one has to ensure that these fluctuations never become strong, which would invalidate 
the linearized approximation. 

A remark is in order. Since we are at finite temperature we have compact euclidean 
time and it is not clear right away what is meant by dynamics and evolution. Our 
coordinates r and x M with /i = 1, 2, 3 are all spatial. To study the "adiabatic" dynamics 
of a single test particle one should add a test-particle time coordinate with g tt = 1, g t/Ji = 
g tr = 0. Indeed, a single test particle only propagates in the thermal background specified 
by the metric but does not interact with the thermodynamical system. Thus in the 
following, we use the five coordinates: test-particle time t and four space coordinates x % 
split into r and x M , /i = 1, 2, 3. 

The relevant part of the wave equation is 

= . (6.1) 

We will soon justify the need to keep only this term in the wave equation for our class of 
examples, where the metric behaves in general as in (|6.4j) below. In order to check the 
normalizability of the two independent solutions of the wave equation, we will use the 
Sobolev norm 



, , v -99 tt ^H + 7; / V^gg'idrfdrf , (6.2) 

2 Js 2 je 

where q is a constant and the integrals are performed on a constant-time hypersurface 
S. The Sobolev norm is bounded from above by the energy of the fluctuation associated 
with the wave, and guarantees that the back-reaction of the fluctuation is indeed small 
37]. According to this criterion there exists a unique self-adjoint extension of the wave- 




operator and a unique evolution of the system, thus rendering the space wave-regular, if 
only one of the two independent solutions of ( |6.1| ) is normalizable near the singularity. 
This solution is then kept, whereas the other is discarded. We should note that this 
criterion is non-trivial for time-like singularities. For null singularities the space-time 
is globally hyperbolic and there should be a unique self-adjoint extension of the wave- 
operator since the evolution is ordinary Cauchy, and hence unique in this case. 



10 This has been shown in the context of the AdS/CFT correspondence and in warped factor com- 
pactifications in |38| (for an exhaustive exposition see |39[ The proof in the present case can also 

be given along similar lines. 
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The choice of Sobolev norm is compatible in our cases with an alternative criterion, 
namely that in geodesically incomplete space-times there should be no leak of conserved 
quantities through the singularity. For the translational isometries that remain unbroken 
by our metric ansatz ( 3.17|) , this amounts to the condition (see, for instance, [ II] ) 



^g rr diijd r ij = (6.3) 

at the singularity. In order to have a unique evolution of the system, only one of the two 
independent solutions of eq. ( |6.1| ) should satisfy (|6.3|). For the other solution, which is 
discarded, the right hand side of ( |6.3j ) should diverge. 

Returning to our metrics, we have seen that they develop a curvature singularity as 
r — > + or r — > r^. We will discuss the case of a singularity at r = + , the other case 
being analogous. Then the asymptotic behaviour of the metric is of the form 

ds 2 ~ dr 2 + (const.) r 2v r\ ilbp dx il dx p , as r — > + , (6.4) 

with v being a positive constant that can be read from the appropriate expressions 



in sections 4 and 5. In the case of the strongly coupled solution ( 5.18 ) around the 
critical point = (0,oOq), we have that v — 1, and for all other cases < v < 1. 

Applying the general formalism for a background behaving as in (|6.4|) we first see that 
the two independent solutions of ( |6.1|) are <~ 1 and ip2 ~ r 1_3z/ (or ip 2 ~ hir if 
v = 1/3). Note that if we had treated the full wave equation in (|6.1| ), we would have 
added a term proportional to r v ip on the left-hand side. Such term perturbs the two 
independent solutions that we have just mentioned by a term proportional to r 2 ~ 2v and 
r 3 ~ 5l/ , respectively. Therefore, near the singularity occuring at r = 0, this extra term can 
be neglected if v < 1, or does not change the power of r if v = 1. 

Applying either one of the above criteria, we find that backgrounds with a single 
singularity behaving as in (|6.4j ) are wave-regular for v > 1/3.0 Hence, for the genuine 
type-II, as well as for the hybrid type-II solutions, the backgrounds are wave-regular 
when the constants c and b, respectively, are negative or zero, while positive c or b are 
unphysical. 

For the heterotic case, we have to distinguish solutions with semi-infinite intervals of 
r (those not going into region 2) and solutions with finite intervals of r (those going into 
region 2). For the former there is only one singularity at r = + and we find that all 
solutions that contain in some limit the region 3 do not satisfy the criteria, and hence 
are unphysical. The only physical solutions in this class are those corresponding to the 
trajectories that connect region 1 to the weak coupling region, or the special trajectory 
that go from the critical point ujq to region 1 or to the weak coupling region. Note that 
for all these trajectories u remains bounded by Uq. Now consider the solutions that have 
a finite interval of r. They have two singularities, one at r = + and one at If at 



11 This condition also guarantees that the solution -02 ~ r 1-3!y that diverges at the singulariry occuring 
at r = 0, is discarded. Hence, if the wave equation that we are examining corresponds to linearized gravi- 
ton fluctuations, these fluctuations will indeed remain small, thus not invalidating the approximation. 
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each singularity one of the two solutions to the wave equation is discarded, no solution 
will be left, as in general the solution discarded at one singularity will not be the same 
as the one discarded at the other. Such a background can be wave-regular only if one 
of the singularities discards one solution and the other singularity does not impose any 
condition, i.e. we need v > 1/3 at one singularity and v < 1/3 at the other. We have 
v = 1/3, 1/7, 1/7 and 1 for regions 1, 2, 3 and the critical point. We conclude that 
solutions corresponding to trajectories from region 2 to region 1 or to the critical point 
are wave regular. The solution corresponding to the trajectory from region 2 to region 
3 however is not wave-regular. It is evident that a better understanding of the string 
rather than the particle propagation is desirable on such backgrounds.^ 



From the string thermodynamic point of view it has been argued [|15| 16, 17] that 
the exponential growth of the number of states forces one to work in the microcanonical 
ensemble, in which case the specific heat may turn negative unless all spatial dimensions 



are compact; of course some compact dimensions may well be large (see also [43 



Within our domain wall ansatz, compactness or not of the space coordinates is not 
an issue since the solutions do not depend on those directions. However, compactness of 
the r coordinate is possible only if the solutions 4>(r), x{ r ) an d u)(r) at some r equal their 
values at some other r. As we discussed in the previous section, this is not the case for the 
type-II solutions. For the heterotic case we have some solutions that are parametrized 
by a finite interval in r. However, without the explicit solution at our disposal, it is not 
clear whether they can be continued to periodic solutions or not. In any case, before 
jumping into conclusions about our solutions being physical or unphysical within string 
theory, one should bear in mind that the compactness criterion was established within the 
microcanonical treatment of string theory, which assumes that the theory remains weakly 
coupled everywhere.^ It is clear from our solutions that they always contain regions of 
arbitrarily strong coupling, thus invalidating the assumptions that could sentence some 
or all of them unphysical. 

Another argument in favour of our solutions is their supersymmetry. As shown above, 
all solutions of the first order equations preserve half of the supersymmetries, by their 
construction, while the standard high-temperature phase of strings breaks all supersym- 
metries. Clearly, at this point we cannot say anything further as definite about the 
relevance and physical properties of our different domain wall solutions within string 
thermodynamics. Ideally, a string inspired criterion would render physical all or a proper 
subset of those backgrounds, which have already been named physical using the criteria 
of supergravity in this section. 



12 The above criteria have been succesfully applied within the AdS/CFT correspondence in many 
backgrounds related to the Coulomb branch of J\f = 4 and J\f = 2 super Yang-Mills theories |Q . In those 
cases the lack of information on string theory beyond the supergravity approximation is compensated 
by the information related to gauge theory expectations. 

13 It is somewhat ironical that the only solutions that perhaps could be continued periodically actually 
are never weakly coupled. 
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7 Conclusions and discussion 



In this paper we have constructed domain wall solutions of the effective supergravity 
which describes all possible high-temperature phases of the known N = 4 superstrings. 
We have used a universal thermal potential that contains the s, t, u moduli and their 
couplings to the three (lightest) winding modes, which can become tachyonic above the 
Hagedorn temperature, thus triggering the thermal instabilities. Our solutions contain 
non-trivial winding fields, which vary with respect to the domain wall variable r, and 
exhibit a common property in that they always extend to regions of strong coupling. We 
have presented the exact solutions for the type-IIA and type-IIB sectors of the theory 
and a certain (self-dual) hybrid sector, and investigated the general structure of the 
solutions for the heterotic sector by extracting their behaviour in the vicinity of the 
weak and the strong coupling points. From a geometrical point of view we found that 
none of the type-II solutions has support on a compact spatial dimension, namely the 
direction parametrized by the variable r. For the heterotic case, the same is true for 
all solutions that contain a weakly coupled region, while solutions that never are weakly 
coupled have support only on finite intervals in the variable r. Whether or not these latter 
solutions can be continued into periodic functions of r could not be answered without 
having an explicit solution at our disposal. From a supergravity point of view, for each 
type of theory, there are subclasses of domain wall backgrounds which are selected by 
imposing specific boundary conditions that lead to consistent propagation of a quantum 
test particle. 

The main property of the domain wall solutions is their supersymmetry; they all 
satisfy first order differential equations, which arise as BPS conditions for gravitational 
backgrounds that preserve 1/2 of the supersymmetries. The supersymmetry of our solu- 
tions is the key to understand why, for all of them, tachyonic instabilities never occur, 
even though the temperature may become arbitrarily high. Of course, we are only con- 
sidering solutions of an effective supergravity and not of the full string theory, but it may 
well be that they point to a new finite-temperature phase of superstrings which is su- 
persymmetric and has no thermal instabilities, i.e., no Hagedorn temperature. One may 
then speculate that, as superstrings are heated up from zero temperature, they prefer to 
go into this more symmetric phase which is stable due to supersymmetry, and that the 
"ordinary" high-temperature phase with the Hagedorn instability is never reached. 

Let us now address a few points one might object to our solutions. First, one might 
worry about the status of the domain wall ansatz in theories of gravity coupled to scalar 
fields and having a potential with run-away directions that account for the usual insta- 
bilities at the Hagedorn temperature. Note however, that this is not a priori forbidden 
by any general considerations alone. Recall that domain wall solutions with infinite 
tension (corresponding to infinite central charge of the supersymmetry algebra) have 
already found numerous applications in supersymmetric gauge theories, most notably 
in supersymmetric quantum chromodynamics with one massless flavor, where one has 
models with run- away vacua. In fact, these are models with rigid supersymmetry which 
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admit stable field configurations that restore one half of supersymmetry, as opposed to 
the stable but non-supersymmetric ground state, and which are characterized by con- 
stant positive energy density [E5 .The domain wall solutions that we have constructed in 



supergravity for the problem of string thermodynamics could be viewed as the gravita- 
tional analogue of such supersymmetric configurations, having their own characteristic 
properties. Their stability is guaranteed by the saturation of the BPS bound and can be 
further supported by analysing the spectrum of the graviton fluctuations on the domain 
wall backgrounds. Based on general principles, one can show that the graviton spectrum 
is obtained by computing the energy levels of an equivalent non-relativistic problem in 
supersymmetric quantum mechanics with a Schrodinger potential that is determined en- 
tirely by the conformal factor of the metric. The details are not important as the only 
relevant point here is that the spectrum is bounded from below (by zero) by the general 
properties of supersymmetric quantum mechanics, thus rendering a physical spectrum for 
the graviton fluctuations on the supersymmetric domain wall backgrounds that preserve 
three-dimensional Poincare invariance. 

Second, one should also study metric fluctuations taking us away from the class of 
conformally flat metrics of the domain wall backgrounds. We have not investigated this 
issue, but we expect that an analysis along the line that can e.g. be found in |39|, |4(| is 
possible. 



Third, a well-known instability of hot flat space was analysed in |46[ . In this simplest 
example it was found that the classical Jeans instabilities arise as a tachyon in the graviton 
propagator, using small fluctuations about hot flat space, and it was further suggested 
that hot flat space will nucleate black holes. Inevitably, this is the fate of gravitational 
systems due to the attractive nature of gravitation. In the case of a relativistic medium, 
the Jeans instability implies that a thermal ensemble of sufficiently large volume will 
collapse into a black hole. Of course, it remains to be seen how these results generalize 
to domain wall backgrounds of our effective supergravity by performing calculations as 
in ordinary quantum gravity. One might fear that most of our solutions will be afflicted 
by Jeans instabilities because they do not have support in small volumes. However, 
their defining BPS property ensures their stability quantum mechanically as well. Put 
differently, the BPS equation which is interpreted as a "no force" condition, although 
this interpretation is more appropriate for asymptotically flat backgrounds with a definite 
Newtonian limit, ensures that no gravitational collapse will occur. 

Clearly, spherically symmetric black hole type solutions of the effective supergravity 
have to be analysed in the future before addressing the long standing problems of quan- 
tum gravity within string thermodynamics. These general remarks put in perspective 
future attempts to construct other non-trivial solutions of the effective supergravity for 
the s, t, u moduli and the string winding modes beyond the supersymmetric domain wall 
ansatz. An important lesson that we already learned, even from the simplest solutions 
we have constructed here, is that they cannot be entirely confined in regions of weak 
coupling only. 
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It would be interesting to study the most general domain wall solution of the effective 
supergravity that describes all possible phases of N = 4 superstrings at finite temperature 
and includes all six scalar fields. The resulting system of first order non-linear differential 
equations does not seem tractable by analytical methods alone, but a combination of 
numerical analysis and analytical calculations around certain points may well provide a 
reasonable global picture, similar to what we did in the heterotic case. In the appendix, 
below, we summarize for completeness some general remarks about the search for periodic 
orbits in order to appreciate the difficulty to establish analytic criteria for their existence 
in the general case. 

Finally, it will also be interesting to study the higher dimensional interpretation of 
our domain wall configurations, as in other theories of gauged supergravity. This might 
be related to thermodynamics of string theory in D-brane backgrounds as has already 



been discussed in the weak-coupling limit in [fig , f49fl . One may also try to extend the 
formalism to construct domain wall junctions. 
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A Remarks on periodic orbits 



Although we have seen quite explicitly that all our type-II solutions are non-periodic 
in the domain wall coordinate r, and the same probably is true also for the heterotic 
solutions, it is much less clear whether this will still be true for the solutions of the full 



on numerical methods, it would be nice to have at least some analytical tools at one's 
disposal. One possibility is to try to establish the existence or non-existence of periodic 
orbits in the solution space of our non-linear system of first order differential equations for 
the domain wall configurations. Clearly if there are such periodic orbits in the solution 
space (e.g. in the <fi — uj plane in the heterotic case) then the corresponding solution will 
also be periodic in r. Note however that the converse is not necessarily true as would 
be examplified by the heterotic trajectories going into region 2, should it turn out that 
they can be periodically continued in r. Of course, then one might still consider them as 
periodic orbits in the (fi — uj plane, although degenerate ones. Finally, for reasons already 
discussed at length in the introduction and discussion section, it is not clear whether 
periodicity may serve as a criterion within string thermodynamics to distinguish physical 
from unphysical solutions. 

It is well known from the theory of dynamical systems that the existence of periodic 
solutions, i.e., closed loops in the two-dimensional parameter space (0, &), of a first 
order system = P(<p,uj) and Co = Q(<f),uj) is a very important and delicate question 
that often can be established only numerically (unless the general solution is known in 
closed form). The boundary conditions select specific orbits 4>(u) by fixing the integration 
constants, and so by singling out the periodic orbits, if there are any among the solutions, 
corresponds to specifying some physical boundary conditions. 

To appreciate the difficulty in establishing the existence of periodic orbits analytically, 
we recall briefly some basic elements of Poincare's theory (see, for instance, [f|7| for an 
elementary account). There, one has the notion of a limit cycle, i.e. a stable closed curve 
in the parameter space, independent of initial conditions, towards which solutions tend 
in an asymptotic sense, or from which they unwind, as it were, as ( —>■ ±oo. A well 
known theorem states that if a limit cycle exists in a given system, then the existence of 
periodic orbits is guaranteed, even though their explicit construction could be a difficult 
task. On the other hand we have Bendixson's theorem stating that if one considers the 



in any given domain in (<f), uj) which is bounded by a simple curve C, the system will have 
no limit cycles inside that region if J has constant sign in it. Therefore, all we can learn 
on general grounds is for which regions there are no limit cycles; if J changes sign in a 
bounded region (and hence it is bound to vanish somewhere in it), the system will not 
necessarily have limit cycles in it, but there is a good chance that it will. Furthermore, 
there are systems with periodic orbits that have no limit cycles. 



set of six equations ( |3.19| ) - ( |3.24| ). Since for the latter one probably has to rely heavily 



function 
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According to these results, we can already see a difference between the type-II and 
heterotic sectors concerning the possibility to have (or not to have) limit cycles in their 
parameter space. Computing Bendixson's function for these two sectors (using the coor- 
dinate C that decouples the x-dependence from the flows of the fields <f> and a;) one finds 
the result 

J„ = -^=(2e^ + e-^(l + 3sinh 2 cu)), 

^= (2e^ + 1 + 2sinh 2 c) , (A.2) 
—j= (cosh(v / 20) - 3sinh(v / 20)sinh 2 cu) . 



Jhyh — 
Jhet = 



It is obvious that Ju and Jh y b are strictly negative everywhere in the parameter space 
(<f), uj), and so according to Bendixson's theorem we learn that there are no limit cycles 
in these cases. Of course, this by itself does not rule out the existence of periodic orbits, 
but we already know from our explicit solutions that for the pure type-II and hybrid 
type-II sectors there are none. On the other hand, the heterotic sector could support 
limit cycles, since Jhet vanishes along the curve 

coth(v / 20) = 3sinh 2 c; . (A.3) 

This curve has two branches, one in the upper (0, uj) half-plane and the other in the lower, 
which are related to each other by the discrete symmetry id — > —uj. Drawing the curve 
in the upper half-plane we see that it starts asymptotically from id = +00 at = and 
drops monotonically to the value uj = lny^ < uJq, which is approached asymptotically 
as — > +00. Therefore, if one considers bounded domains in the parameter space which 
intersect with Bendixson's curve Jh e t = 0, they will potentially contain limit cycles, 
although their existence is not at all guaranteed. Note that this curve is crossed precisely 
by the trajectories going to region 2. Of course, these trajectories are not contained in 
any bounded domain. Nevertheless, this is still suggestive that maybe these trajectories 
can indeed be periodically continued. 

It will be interesting to look for closed trajectories in the higher dimensional parameter 
space with all the relevant fields turned on. 
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